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The images of many distant galaxies are displaced, distorted and often multiplied by the presence 
of foreground massive galaxies near the line of sight; the foreground galaxies act as gravitational 
lenses. Commonly, the lens equation, which relates the placement and distortion of the images to the 
real source position in the thin-lens scenario, is obtained by extremizing the time of arrivai among ali 
the nuli paths from the source to the observer (Fermat 's principle). We show that the construction of 
envelopes of certain families of nuli surfaces consitutes an alternative variational principle or version 
of Fermat's principle that leads naturally to a lens equation in a generic spacetime with any given 
metric. We illustrate the construction by deriving the lens equation for static asymptotically fiat 
thin lens spacetimes. As an application of the approach, we find the bending angle for moving thin 
lenses in terms of the bending angle for the same defiector at rest. Finally we apply this construction 
to cosmological spacetimes (FRW) by using the fact they are ali conformally related to Minkowski 
space. 

I. INTRODUCTION 

It is the purpose of this note to point out, study and apply an idealized construction of gravitational lens equations 
that is of potential use in many physical situations - from exact lensing to the weak-field thin-lens scenario - by means 
of an alternative version of the standard usage of Fermat's principle. 

The fundamental aspect of gravitational lensing theory is the construction of the past light-cone of an observer. 
This directly leads to the idea of the mapping from the spacc of images - the celestial sphere of the observer - to 
the space of the sources - usually to the "source piane", though this specialization is by no means necessary. The 
mapping is carried out by following, backwards in time, the nuli geodesics of the lightcone, from the observer to the 
source. In other words, by knowing where an image appears on the observer's celestial sphere and knowing the nuli 
geodesics that generate the past nuli cone, one could in principle follow the rays back to the source. In addition, 
it is often of considerable importance to know the transit times between the emission of light and its arrivai at the 
observer. In fact, in view of Fermat's principle, the actual path taken by light is a locai extremum of the transit time 
of ali possible neighboring nuli paths, which leads to the formulation of gravitational lensing via Fermat's principio. 
Much of contemporary lensing theory is based on the construction, on a simple background (either Minkowski or 
a cosmological spacetime), for weak fields and with a thin-lens and small-angle approximation, of an appropriate 
transit time function. Then, by the locai extremization of the time function, a lens equation is constructed. In the 
usuai approach, the time function (referred to as a Fermat potential) represents the transit time along ali possible 
nuli curves, not necessarily geodesie, that connect the source and the observer The extremization produces a 

suitable lens equation by selecting those curves that are geodesie, Le., the variation of the travcl time with respect to 
the paths, restricted to the condition that the paths bc nuli, is equivalcnt to the geodesie equation with nuli tangents. 

We wish to show that there is an attractive alternative (and signihcantly different) variational principle (an al- 
ternative Fermat's principle) for use in lensing that can be applied, at least in principle (and, with approximations, 
also in practice), in a generic situation. The basic framework is to begin with a general four-dimensional Lorentzian 
spacetime where the geometry (the metric) of the spacetime is to be considered as the "gravitational lens" . In other 
words, the detailed lens properties are to be coded directly into the metric tensor. We then consider a 2-parameter 
family of nuli (or characteristic) surfaces passing through an observer's world line at a given time. In fact, we assume 
that the family of surfaces is sufficiently generic for the nuli normals at the observer to be distinct and span the 
sphere of nuli directions (often just an open neighborhood of the sphere is sufficient). This family of nuli surfaces 
then contains ali the points on the observer's lightcone (or the open neighborhood), since at the intersection of the 
observer's worldline with each surface the normal vector to the surface is nuli, geodesie, and lies on the surface. Each 
nuli geodesie that passes through the observer's world line on each of the nuli surfaces can then be followed into the 
past. These are the rays that an observer sees and constitute his celestial sphere (or an open neighborhood of it). 



1 



Now consider a point source of light moving along some given (source) timelike world line. We are interested in 
those nuli geodesics moving back in time from the observer to the source, Le., those geodesics that the observer "sees" 
as coming from the source. At any given observer moment, the observer will "see" a number of different light rays 
(or images), each of which, in general, will have intersected the source (or equivalently, will have been emitted by 
the source) at different source times. These countably few nuli geodesics lie on countably few of the nuli surfaces in 
the family that we are considering. Thus there are a number of surfaces in the 2-parameter family that intersect the 
worldline of the source at a point that can be connected to the observer by a nuli geodesie on the surface. However, 
at any observer moment, ali (or almost ali) the other surfaces in the family intersect the source's world line as well, 
at varying times. In general, however, there are no curves lying on these surfaces connecting the source intersection 
point to the observer, that are nuli geodesics. In fact, most of the curves that can be used to connect the source and 
the observer on each of these nuli surfaces are pieeewise spacelike and/or nuli. For this reason, we prefer to drop any 
reference to paths, and keep the argument in terms of the nuli surfaces. The source time t at the intersection point is 
thus a function on the sphere of nuli surfaces intersecting the observer's worlddine at a particular observer's time To: 

t = T(6,<p,T ), (1) 

where (9, ip) label the nuli surfaces. We now ask for the locai extremes of T as a function of (9, <p): 

d e T = T, e {9,<p,T ) = 0, (2a) 
d v T = T )V (d,<p,To)=0. (2b) 

This operation picks out those nuli surfaces that possess curves from the observer traveling backwards to the source 
which are nuli geodesics. It constitutes our version of Fermat's Principle. This version differs from the usuai one not 
merely in form. If one thinks of it in terms of paths connecting the source and the observer, this version of Fermat's 
principle allows for curves that are neither nuli nor geodesie. It concentrates on and varies the nuli surfaces rather 
than the nuli curves. 

For this reason we think of our version of Fermat's principle as an alternative to the usuai one. Correspondingly, 
we refer to the function T as a generalized Fermai potential. With respect to the underlying meaning of our version 
of Fermat's principle, one can see that Eqs. (||) are equivalent to the construction of the envelope of the nuli surfaces 
passing through the observer, which in turn is the past lighteone of the observer. Thus we arrive at the observer's 
lighteone starting from surfaces in a way analogous to the usuai approach, which arrives at the observer's lighteone 
starting from nuli paths. 

In Section II we describe this construction in greater detail and justify the claim that it does pick out the nuli 
surfaces so that nuli geodesics connect the observer with the source thus yielding the past lighteone. In Sec. Ili we 
illustrate the construction for the trivial case of Minkowski space without a lens, while in Sec. IV we illustrate it for a 
Minkowski space background with a static thin lens, the conventional scenario. In Sec. V we apply the construction 
to thin lenses that are moving in order to obtain the corrections to the lens equation due to the motion of the lens. 
In Section VI these ideas are applied to the FRW cosmologies using the fact that they ali are conformally related 
to Minkowski space. We conclude, in Section VII, with remarks and an outlook of the potential applications of the 
construction. 

We find this method quite enlightening concerning the underlying ideas in gravitational lcnsing theory. And 
although it might be difucult to apply in many practical physical situations without the use of heavy approximations, 
preliminary calculations suggest that, with further development, it could well be of practical use. 

II. IMPLICITLY DEFINED GENERALIZED FERMAT POTENTIALS 

We assume that we have a four-dimcnsional Lorentzian manifold, (M, x a ) with a given Lorentzian metric, g a b- 
We further assume that nuli geodesics and nuli surfaces (solutions to the eikonal equation) can be constructed, in 
some locai coordinates, wherever needed. In general, nuli surfaces develop wavefront singularities and, in principle, 
one must worry about them. However, the trouble caused by the singularities is easily bypassed if the surfaces are 
parametrized by the nuli geodesics that rule them. The fundamental objects for us will be complete integrals of the 
eikonal equation: two-parameter families of foliations of the spacetime by nuli surfaces, so that at any spacetime point 
there is a sphere's worth of nuli surfaces through that point. The two parameters are arbitrary, but for our purposes 
we specify them as points on the sphere. In the following, then, we adopt this particular choice of parameters to 
label the solutions of the eikonal equation in the complete integrai, and we specify them by complex stereographic 
coordinates (£, £), rather than the regular spherical coordinates (9, (p). The complete integrai is given analytically by 
an expression of the form 
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Z(x a ,CO (3) 

such that for each fixed value of (£, £) the level surfaces of (|^) are nuli. Equivalently, (Q) satisfies the eikonal equation 

g ab d a Zd b Z = 0, (4) 

with two further conditions: the four functions (Z,d^Z,d^-Z,d^Z) forni a rank-four set with respect to x a almost 
everywhere, whereas globally, the three functions (Z,d^Z,d^Z) forni a rank-three set. How these conditions are used 
is explained shortly below. In addition, the nuli covectors Z a (x a ,£, () = d a Z(x a ,(,C) a t fixed x a sweep out ali nuli 
directions as (£, vary, thus reproducing the locai lightcone. We rewrite (^) as the two-point function 

G(x a , x a , C, C) = Z(x a , C, - Z(x a 0: C, C) - 0. (5) 

Any two points (x a ,xg) satisfying Eq. (|5|) for fixed (CO ue on t ne sanie nuli surface. It is clear that the gradient of 
G with respect to either point is a nuli covector. If xg is fixed, however, then the points x a that satisfy Eq. (||) for 
each (C lie on different nuli surfaces, but xg belongs to ali. So we think of xg as a special point where ali the nuli 
surfaces intersect. 

The point xg is chosen to lie on the observer's worldline given by xg = xg(To), Tq being some observer time. Though 
it is not needed and is used only for convenience, we choose, in the neighborhood of the source, a special coordinate 
system where the source remains at a fixed spatial point, x l , and only the source time x° = t evolves. The equation 
G(x a , Xq, () — then becomes 

G(x i ,t,xZ(T ),(,()=0 (6) 
or, rewritten (since dtG ^ by assumption) as 

t = T(x i ,T ,(,Q, (7) 

it gives rise to our time function, Eq. ([!]). As a matter of notation, by analogy with the usuai theory of gravitational 
lcnsing it is naturai for us to refer to T as a generalized Fermai potential. However, the function G arises more 
naturally from the starting point of the eikonal equation, therefore, with a slight abuse of terminology, we refer to it 
as the ìmplìcìt Fermai 'potential, and for ali practical purposes we use Eq. (^|) instead of Eq. (|7|) . 

The extremal values of T are calculated in terms of the implicit Fermat potential G by implicit differentiation of 
Eq. (Q). Thus d^T and d^T are obtained from 

G, t d c T + d c G = 0, (8a) 
G, t d- c T + <9 f G = 0, (8b) 

so that the vanishing of d^T and d^T implies the vanishing of both d^G and d^G, and the reverse, i.e., 

(d c T, <9 C -T) = O (d c G, 9 C -G) = 0. (9) 

Thus, setting OqG = d^G = gives us the extreme values of T. There is a deeper meaning to the extremization of T, 
in connection with the nuli surfaces in the family, as we subsequently show. 

Returning to Eq.(||), with fixed xg, we see that we have a two-parameter family of surfaces through a fixed point xg. 
A new surface, with common tangent curves ali passing through xg, is the envelope of the family and is constructed 
by requiring that d^G = d^G = 0. The triple 

G(x a ,x§,C,C) -0, (IOa) 
d c G = 0, (lOb) 
<9 C -G = 0, (lOc) 

defines the envelope. For x a ^ xg and in a region without wavefront singularities f|,|]], Eqs. (lOb) and (lOc) can be 
algebraically inverted so that 

C-T(x a ,xg), (Ila) 

C-T(x a ,x°). (llb) 

When they are substituted into (|l0|), we have the equation for a three-surface, namely the envelope: 
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G e nv{x a i Xq) = G(x a , Xq, T(x°, Xq), T(x°, Xq)) — 0. 



(12) 



That G en y ( x a , X n) = is a nuli surface follows from both the fact that dG/dx a is a nuli vector (by assumption) and 
from ( |l0b| )-( [l0c| ). That it is the nuli cone of the point Xq is shown by demonstrating that at x a = x" the surface has 
a sphere's worth of tangents. This is seen by the following argument: in the limit as x a 



identically satisfied for ali values of (£, Q. This is a degenerate point where Eqs.flllD do noi hold. It then follows that 
the tangent vector, given by 



Eqs. (IlObD-fllOq 1 ) are 



dG e 



dx a 



(13) 



has multiple values, yielding a nuli vector that sweeps out the tangent-space nuli cone at Xq, as (C,C) nioves over the 
sphere. 

Our task is now to show that this extremization process picks out the particular surfaces that connect the source 
and obseryer by nuli geodesics; namely, the values of (£, Q for each pair (x a , Xq) given explicitly by ( |TT| ) - or implicitly 
by ( |l0bj )-( [l0^) - label a geo d esie nuli vector whose integrai line passes through Xq and x a . In order to show this, we 
notice that when Eqs. (|l0b|)-fl0c|) can not be inverted, Le., when the Jacobian 



d C( G d ?r G 



d^G %G 



(14) 



vanishes, the nuli surfaces develop wavefront singularities M. But by the assumption that Z(xq, x a , (, () is a complete 
integrai and the rank conditions, in either case the three equations ( |l0| ) can be solved for three (say x l ) of the four 
spacetime coordinates x a , in terms of a fourth one (a:*) and (C)C)j namely: 



x* =X i (x%,x*,t,0- 



(15) 



As a consequence of (|To|) , it turns out that the curves described by ( |l5| ) by keeping (Ci C) hxed are nuli geodesics. 
This can be seen by implicit differentiation of (|l0| ) with respect to x*. Defming t a = (l,dX l /dx*\<* j), we obtain 



G, a (x a ,x a ,t,0t a =0, 
d c G, a t a = 0, 
d?G, a t a = 



(16a) 
(16b) 
(16c) 



where {G, a , d^G, a , d^G ia ) are independent, by the rank condition, again. Eq. (16a) implies that t a has no component 
pointing out of the surface, so t a — ag ab (G,b +Pd^G,b +(3d^G,b)- Using this in Eqs. (16b) and ( |16cj ) implies that 
j3 — (3 — and hence t a — ag ab G,b, and thusthe tangent vector is proportional to g ab G,b being therefore nuli, and the 
integrai curves are nuli geodesics. 

We thus obtain the parametric description of ali the nuli geodesics through Xq . Notice that having this parametric 
description of the lighteone (in terms of the observer's celestial sphere) is entirely equivalent to having a lens equation 
relating the angular position of the source at a given distance to the angular position of the image on the celestial 
sphere, plus a time of arrivai equation yielding the transit time of the light signal from the source to the observer 
||[7]]. To see this, the point Xq is taken to represent a point on the observer's worldlinc. By treating (C>C) as the 
celestial sphere of the observer, x* as a measure of "distance" backwards along the nuli geodesie to a source and x l 
as representing the "time" and "angular position" of the source at the moment of emission, we have constructed a 
lens equation together with the transit time equation. Thus one can see that, in a completely general context, our 
version of Fermat's Principle, Eqs. (|2|), leads immediately to the construction of the observer's past lighteone, with 
associated lens and time of arrivai equations. 

When we extremize the function T we could obtain a maximum, a minimum or a saddle point. We wish to know 
how to rephrase the conditions for maximum, minimum or saddle in terms of our implicit Fermat potential, G. In 
principle, we need to look at the eigenvalues of the matrix containing the second derivatives of T evaluated at the 
extrema. We now switch to real variables (u, v) instead of (£, (), via ( = u + iv , and calcitiate the eigenvalues of the 
matrix 



S 2 T: 



T T 



T T 



(17) 
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Because of T being implicitly defined via (g), we just need to take second derivatives of (g). For instance, if we take 
two u— derivatives of (g) (with (0) for t) we obtain 

G,uu +T, U (2G, u t +G,tt T, u ) 4- G,t T, uu = 0. (18) 

Sirice we are interested only in the value of T, uu at the extremum, then T, u — and 

T, uu =-^, (19) 

which presents no difficulty, since by assumption G,t must be nonvanishing. Likewise, the other components of the 
matrix come out proportional to the second derivatives of G: 

GyuV \ 

(20) 

Thus, up to a factor of — l/G,* the eigenvalues are the same as those of S 2 G. In order to fix unnecessary sign freedom 
(since — G for our purposes is just as good as G), we assume that G,t < at the extremum in question. It remains to 
rcphrase our conditions for maximum, minimum and saddle in terms of the complex variable Since £ = u + iv, we 
have d/du — d/d( + d/d( and d/dv — ì(d/dQ— d/d(), thus 

G, U u = G,^q +G,^ -\-2G,qq , (21a) 
G, vv = — G,£f — G,££ +2G,^ , (21b) 
G, uv =i(G, Cf -G, f - f ). (21c) 

If we think of G as spin-weight zero function on the sphere then, using the envelope condition, G,^ = G,^ = we have 

5 2 G = (1 + CC) 2 G, CC , (22a) 

>cc- 

and hence 

9 / 99G + iìe(g 2 G) -Im{d 2 G) 

8 2 T=— 1- |. (23) 

(1 + CC) G, t y _/ m (g2 G ) g g G _ i?e (g2 G ) 

The eigenvalue equation for S 2 T is given by 

(À - 9§G) 2 - 9 2 G§ 2 G = 0, (24) 

where A = — (1 + (() 2 G,t A/2 and A is an eigenvalue of S 2 T. Since G,t is assumed to be negative, then the sign of A 
is the same as the sign of A. The solutions are 

A± = 9§G± |3 2 G|. (25) 

Finally, the conditions are: 

• if |35G| < |8 2 G|, the extremum is a saddle, since the eigenvalues have opposite signs. 



5 2 G= (1 + CC) 2 G, C - C -. (22b) 



if 39G > |9 G|, the extremum is a minimum, since both eigenvalues are positive. 



• if S3G < — |9 2 G|, the extremum is a maximum, since both eigenvalues are negative. 

On the other hand, the envelope develops singularities when 

mG 2 = S 2 G§ 2 G (26) 

as anticipated earlier (the vanishing of the determinant of the Jacobian matrix) . The type of extremum is of funda- 
mental importance. For instance, in a very weak field, no singularities develop in the observer's lightcone, there is 
no multiplicity of images and the lens equation represents a locai minimum of the Fermat potential. Increasing the 
strength of the field, in a situation where there are three images, one will be a minimum, one a maximum and the 
remaining one will be a saddle. The minimum corresponds to the primary image, produced by a ray that does not 
encounter a caustic in its path. The saddle corresponds to a ray that passes the caustic once. The maximum yields 
the ray that passes the caustic twice and yields the faintest image. The generalization of this interpretation to cases 
of more than three images is complicated and lies beyond our present interest. 
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A. Parametric Version of Implicit Fermat Potential 



We finish this section by presenting a very important mcthod for describing implicit Fermat potentials. In many 
cases it is difScult or impossible (as we will see) to obtain a closed-form G for a family of surfaces G(x a , xg, C, C) = 0> 
especially if the surfaces self-intersect and have singularities. In such cases, it is simpler to describe the family of 
surfaces in parametric form adapted to the nuli geodesics ruling the surface, namely, by specifying the map 

x a =r a (s,r,q,<:,0 (27) 

where (s, r, q) = s l are coordinates on the surface, given for each fixed value of (£, £). Since each member of the family 
(for fixed (£, ()) is a nuli surface, it is automatically ruled by nuli geodesics (except at singular points of the surface). 
Therefore the surface coordinates s l can always be chosen as geodesie coordinates: (r, q) labeling the nuli geodesics 
ruling the surface, and s as an affine parameter along the geodesics. However, our interest is in surfaces that intersect 
at a given point x$. Therefore, ali the surfaces in the family must contain the common point Xq. We can make use 
of this fact by adapting the coordinates (s, r, q) to this point. Since s is a parameter along the geodesics, we fix it by 
demanding that the 2-surface s = be transverse (to the nuli geodesics in the nuli surface) and contain the point x$. 
This is, then, the initial-data surface for the geodesics in the nuli surface. On this initial-data 2-surface, the origin of 
coordinates for r and q can be taken to be x$. In this manner, ali the nuli surfaces in the family will single out the 
observer's point as the point with s = r = q = 0, namely: 

4 = r Q (o,o,o,c,C) (28) 

so that r a depends on x$, i.e., 

x a = r a ( S ,r,q,x a ,(,(). (29) 

Oftcn we take x$ as the observation time To at the spatial origin and omit it from the equations. 

In the case that the family of surfaces G(x a , Xq, C, C) =0 is given parametrically via (^7), the envelope construction 
proceeds in an implicit manner as follows. In the first place, the ruling of the nuli surface via nuli geodesics guarantees 
that the Jacobian matrix of the map (|27|), namely dT a /ds l , has rank three (except at singular points). Therefore, 
locally, from the equations ([2?]) one can always chose 3 equations that can be solved for s\ The fourth equation 
remains a function of s l . So equations ( |2^ ) can always be viewed in the form 

r( s \C,C;^)-^ = o, (30a) 
iV,c,C;^)-z 4 = o, (30 b) 

r = r 4 (30c) 

wher e ( |30b| ) is interpreted as the equation G(x a , xg, C) — if s * is thought of as implicitly given in terms of x l via 
(|30a|) . In other words, 



G{x\ xl c, C) = IW 1 , C, C; x a ), c, C; zg) - * 4 - (31) 

In the envelope construction we then have 



= d c G\ xa = d ( T\ s 



dcs%, (32) 

c 



where dQS l \ x i is determined by implicit differentiation of (30a). Both Eqs.(j3l]) and (|32| ) will be extensively used in 
the following sections. 

This appears to be a convenient method for constructing implicit Fermat potentials and for obtaining lens equations, 
as demonstrated in the following sections. 



III. THE IMPLICIT FERMAT POTENTIAL IN MINKOWSKI SPACETIME 



The case of Minkowski spacetime is quite trivial but is helpful as an illustration of the main ideas underlying the 
concept of the implicit Fermat potential G. In cartesian coordinates x a , nuli surfaces in Minkowski spacetime are level 
surfaces of the solutions to the eikonal equation r] ab G, a G,b = 0, with i] a b =diag(l, —1, —1, —1). There is a 2-parameter 
family of distinct foliations by nuli planes with the two parameters corresponding to a sphere's worth of nuli directions 



G 



at any given spacetime point. Using complex stereographic coordinates, the nuli directions can be specified by / a (C, C) 
given by 

' (C ' C) - V2 \ ' 1 + CC' i + CC ' 1 + CC7' 1 j 

The nuli planes with nuli normal l a containing the observer's location Xq(Tq) are given by 

G(x a , Xq, C, = (x a - Zq%(C, = 0. (34) 

The envelope of this 2-parameter family of nuli planes is the lighteone of the point x$ and is obtained by taking the 
partial derivatives with respect to C and C °f the function G(x°, xg, (, and setting these derivatives equal to zero. 
In this case: 

^-(- a --S)^ = o, (sto) 

^ = (* a -*S)^=0. (35b) 

Eqs. ([54]) and ( |35| ) define a 3-surface in the spacetime. One way to obtain a single equation defining such a surface is 
to solve ([35j) for £ as a function of x a , wherever possible, and substitute ((x a ) back into (HJ). Then the envelope is 
thus given, in principle, by 

G env (x a ) = G(x a ,x a Q ,((x a ),((x a )) = (36) 



In this case, this procedure can be carried out in closed form. Using the notation A X — X ' 

Eq. (p5bj) is 

equivalent to 

(Ax - iAy)C 2 - 2Az( - Ax-ìAy = (37) 

with solution 

„, Az±Ar , 

Ax — lAy 

Ar = V(Ax) 2 + (Ay) 2 + (Az) 2 . (38b) 

(The other eq uati on, ( |35a| ) , is the complex conjugate of (|35b|) and yields no new information). Substituting C( xa ) 
given by Eq. (38a) and its complex conjugate C(x a ) into Eq.fl34|) yields 

G env (x a ,x%) = G(x Q ,xg,C(x a ),C(x a )) = ^=(A*T Ar) = 0, (39) 

where the minus sign is for the future lighteone and the plus sign for the past lighteone of the point Xg . Consequently, 
the set of three equations (37) and ( p4] ) are equivalent to the lens and time of arrivai equations. Because there is no 
timo occurrence in (|3"7f), then one can think of it as the lens equation since it relates the spatial position of source 
(x, y, z) to its corresponding image on the celestial sphere (CjC)- I n principle, solving for two of (x,y, z) in terms of 
(£, C) and substituting into Eq. ( |34| ) turns Eq. ( |34| ) into an equation relating the time of emission t from a source at 
(x, y, z) with image direction (£, £) seen at the observer xg. 



IV. THE IMPLICIT FERMAT POTENTIAL FOR A STATIC ASYMPTOTICALLY FLAT THIN-LENS 

SPACETIME 

In geometrical terms, a thin lens can be modeled as a spacetime with vanishing curvature everywhere except 
at a timelike surface, where the "lens" lives. Accordingly, the spacetime consists of two fiat spacetimes, matched 
appropriately at the lens' worldshcet, with the observer lying in one fiat spacetime and the source lying in the other 
one. We refer to these two fiat spacetimes as the observer's side and the source's side of the spacetime. 

It is clear that the familics of nuli surfaces to be used on the observer's side of the spacetime are ali the nuli planes 
that contain the observer's location, labeled by the direction of the nuli normal. On the source's side, however, the 
appropriate nuli surfaces to use are not nuli planes but the distortion of the nuli planes by the lens. For this reason 
they must, in general, be given in parametric form. We use the parametric method, via Eq. (|27j), explained in Section 
II, for the envelope construction. 

For the purpose of clarity of presentation, we first develop our construction of the nuli surfaces in a lower-dimensional 
setting. We then generalize the argument to 3+1 dimensions. 
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A. 2+1 static lens planes 



The 2 + 1 lensing setting is describcd in cartesian coordinates x a = (t, x, y). The observer is at = (To, 0, 0), and 
the lens surface is a piane located at y = yi. The metric of the spacetime is fiat everywhere except at the lens piane: 

ds 2 = c 2 dt 2 - dx 2 - dy 2 . (40) 

Our aim is to construct a 1-parameter family of nuli foliations which in principle can be given as the level surfaces 
of a function G{x a ,(p) which replaces the G(x a ,(iO 01 the four dimensionai lensing scenario. (The origin, xg, is 
omitted.) This function will be entirely equivalent to the Fermat potential in this lensing scenario, and the envelope 
of the surfaces G(x a ,ip) = will simultaneously give us the lens equation and time of arrivai. In the approximation 
of small angles, the Fermat potential and the lens equation will reduce to the standard astrophysical lensing scenario. 

Roughly our construction is as follows; we define a 1-parameter family of nuli surfaces, (parametrized by ip) by 
imagining a straight line of photons traveling (backwards in time) parallcl to each other (nuli rays) in a direction 
making an angle tp with the y— axis. As the photons move, they trace a nuli piane in the spacetime. The moment 
each photon arrives at the lens piane, it is "detained" for a length of time T that depends only on the point at which 
the photon hits the lens piane. After this time, the photon is released in a new direction tp' = ip — a, deviating by an 
angle a from its originai direction. The bending angle a depends only on the point at which the photon hits the lens 
piane, and it should be entirely determined by the "detention time" T . Because different photons arrive at the lens 
piane at different times and are detained for different lengths of time, the wavefront they make up after leaving the 
lens piane is not necessarily piane. 

The scheme that we have just described has a corresponding analog in the wavefront distortion produced by a 
timi shcct of glass with a variable index of refraction. A lightray impinging nearly normally on the shcet will bencl 
inside the glass due to the component of the gradicnt of the index of refraction perpendicular to its path, exiting 
the glass with a finite deviation angle. A piane incident wavefront then comes out of the glass distorted. The locai 
index of refraction at the point of incidence multiplied by the thickness of the glass is a measure of the time spent 
by each lightray inside the glass, or equivalently, the "detention" time of light inside the glass. For later purposes, 
it is important that we make following point with respect to the limit in which the thickness of the glass tends to 
zero. The "detention time" in the region occupied by the glass defined in this way has two contributions: the transit 
time in vacuum through the region (the "vacuum" time: c _1 x (thickness of the glass)) plus the amount of time that 
light is delayed by the slowing-down in the presence of the medium, with respect to vacuum. In the limit in which 
the thickness of the glass tends to zero, the first (vacuum) term vanishes and the "detention time" coincides with the 
time delay with respect to vacuum. 

In the case of interest to us, the detention time T is determined by the gravitational field of a deflector on the lens 
piane and it coincides with the gravitational time delay with respect to vacuum: 

T = ~J J Udl ( 41 ) 

where U is the Newtonian potential of an isolated mass distribution and the integration takes place along the nuli 
ray or line of sight. The situation we are contemplating is the limiting case in which the support of U tends to zero, 
while the integrai of U remains Constant, namely: U is a distribution along each nuli ray. 

More precisely, on the observer side of the spacetime we consider the parallel lightrays traced by a piane wavefront, 
moving (backwards in time) at an angle ip with the optical axis (the y-axis) as defining our nuli surface. See Fig. [|. 
The wavefront passes by the observer at time To. The (perpendicular) distance from the observer to the wavefront 
at later times divided by c, which is denoted by r, can be used as a parameter along each lightray. The individuai 
lightrays at r = 0, can be labeled by their distance r from the observer. In terms of these parameters, the time 
coordinate of the lightrays is given by 

t = T -r, (42) 

whereas the space coordinates of the lightrays are given by 

x = ct sin ip + r cos tp, (43a) 
y = ct cos ip — r sin ip. (43b) 

(Note that as r increases time is going backwards.) It is clear that r = r = defines the observer at (x, y) = (0,0). 
Each lightray (for fixed (p) reaches the lens piane at r = r; when y = yi, hitting a lens point x — xf. 
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X\ = cri sm tp + r cos tp, 
yi = cri cos tp — r sin tp. 



These, in turn, can be solved for 77 and r and yield 



cti = xi sin ip + yi cos tp, 
r = xi cos tp — yi sin ip. 



(44a) 
(44b) 



(45a) 
(45b) 



Eqs. @ hold while r < 77. 

On the other side of the lens line, the wavefronts will be determined by two things: the bending angle a(xi) and 
the gravitational time delay imposed on the lightrays at the lens line, which we here denote by 7~(xj). The bending 
angle is to be determined in terms of T(x;). The bending angle and the time delay can be incorporated into the 
parametric equations for these lightrays in the following way: 



x = c(r — ti — T) sìn{tp — a) + xi, 
y = c(r - ti - T) cos(<yS - a) + y h 



(46a) 
(46b) 



These equations hold for r > ti + T. For times in between, the lightrays are "delayed" at the lens piane, so we have 



x = Xi, 

y = vi, 



(47a) 
(47b) 



for ti < t < ti + T . Eqs.(^3[), (pq), and ( ]47|) represent the parametric expression of a one-parameter family of nuli 
surfaces which in principle can also be expressed by a single equation of the form G[x, y, t, tp) = away from caustics. 
On the observer side this can be done explicitly by eliminating r from Eqs. ( [Ì3|) yielding 



G(x, y, T,tp) = y cos tp + x tan ip — r = 0. 



(48) 



On the source side, however, this can not be done explicitly. The parameter that labels the nuli geodesics is xi, 
as op pose d to r on the observer side. They are related by Eq. ( [45b| ) . The idea would be to eliminate xi between 
Eqs. (46a) and ( 46b ) but t his is impossible since both T and a are functions of x/. We thus need to define xi as a 
function of x implicitly by (46a). We take (46b) as the equation that defines the function G where xi is the implicit 
function of x given by (46a): 



G(x, y,T,tp) = y — yi — (ct - xi sintp - yi cos ip - cT(x;)) cos(<p - a(xi)) = 0. 



(49) 



We have used (45a) to eliminate r; in terms of xi and yi. 

The envelope is obtained by setting dG(x,y,t, ip)/dip — 0, where the partial derivative is taken at Constant values 
of (x, y, t). On the observer side, Le., for y < yi, taking d/dip on (E§) we obtain 



y sin ip + x cos tp = 



(50) 



x 

tanip = —, (51) 

y 

an equation for tp as a function of x a on the observer's side. Substituting ( |5Ì| ) back into ( ^8| ) we have (for y < y{), 

(52) 



^x 2 + y 2 + c(t - T ) = 

which is the observer's past lightcone, as expected. Notice that when ( |5l| ) is evaluated at the lens piane it yields 



taniy? 



Xl 

Vi ' 



(53) 



From Eq. (|45b|), we see that this implies that r = 0, Le., that these rays come from the observer. In other words, of ali 
the lightrays at angle tp, the envelope condition picks the nuli ray that passes through the observer. This lightray hits 
the lens piane at x = x/ and subsequently bends through an as yet undetermined angle a(xi). In order to determine 
the bending angle a we use the envelope condition on the source's side, imposing on it the condition (|53|). 
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The envelope is obtained by setting to zero the implicit differentiation of G(x,y,T, tp) from Eq. ( p9| ) with respect 
to tp keeping x a fixed. We obtain: 



= sin(yj — a) 



dxi da 
dtp dxi 



1 



(cXq — et — xi sin tp — yi cos tp — cT) 



+ cos(p — a) 



yi sin ip — xi cos ip 



dxi ( dT 



dtp \ dxi 



sin tp 



The quantity |p- is calculated by implicit differentiation of (46a) with respect to tp, keeping x a fixed: 



dxi 
dip 



A 

B' 



(54) 



(55) 



with 



A = cos{(p — a)(cT{) — et — xi sin ip — yi cos tp — cT) 
+ sin(tp — a) (yi sin tp — xi cos tp) 



da 



B = 1 — cos( tp — a)(cTo — et — x\ sinyj — y\ cos tp — cT)— — 

dxi 

^ ( dT 
— sin tp — a) c— h sin 02 



(56a) 



(56b) 



From Eq.(p5|), and some algebra, Eq. ( p4[ ) becomes 
= 

+(yi sin^s - x ; cos tp) 



, 97" 

sm(iy9 — a) + c— h sinw 

da;; 



cos(t^ — a) 



9a 



(cTo — et — xi sin tp — yi cos (/? — cT) da;; 
Using Eq. (|53|), the first parenthesis in the last term of Eq. ( |57| ) vanishes, so that Eq. ( |57| ) reduces to 

5T 



= sin ^9 — sm(tp — a) + e 



dxi ' 



(57) 



(58) 



This is the basic relationship between the bending angle and the gravitational time delay that holds for large as well 
as small angles. Using the regime of small angles, Eq. (F)8|) reduces to 



dT 

a = -e—, 
dxi 

the standard relationship in the astrophysical approach to lensing in a two-dimensional scenario. 
The standard lens equation is obtained in the following manner. Elimination of r — ti — T in Eqs. 



(59) 



) implies 

(x - xi)cos(tp - a) - (y - yi) sin(tp — a) = 0. (60) 
With small angles, and with yi = Di and y = D s (the source distance), this becomes 

x - xi - (D s - Di){tp - a) = 0. (61) 
From the envelope condition, Eq. (|5"3|), with the small-angle approximation (tp = tan^ = xi/Di), we finally have 

x xi (D s - Di) 



D s 



D t 



D s 



-a. 



(62) 



the standard lens equation, which is understood as a map from xi on the lens line to x on the source line (i.e., xi — > x), 
with fixed values of Di and D s . 

As an illustration, Figure^ shows how the envelope construction works for the case of a bending angle that does 
not depend on the point on the lens piane: a(xi) = const. and cT = —axi. The top panel shows a spacetime picture 
of one member of the 1-parameter family of nuli surfaces that is used to construct the envelope. In this case, two 
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nuli planes in the directions tp and (p — a are matched at the lens' worldsheet. The vertical separation between the 
intersections of the two nuli planes with the lens' worldsheet represcnts the time delay T. A spatial projection of the 
resulting nuli surface shows the bending of lightrays, whereas the time delay is observed in the shift of the wavefronts 
in passing through the lens line. The envelope of the ip— dependent family of such nuli surfaces yields the lightcone 
of the observer on the observer's side of the lens, which is matched to another lightcone with a shifted origin on the 
othcr side of the lens. 

A more physical example is provided by the use of a time delay of the form T(xi) = 7ó(l — ln(l + (xi/xo) 2 )). This 
reproduces qualitatively the time delay of a compact lens in the vicinity of the optical axis, since the delay is greater 
for smaller impact parameters. The spatial projection of the nuli surface obtained in this case is shown in the top 
panel of Fig. ||. The center panel of the figure shows the projection of the envelope. In this case, a caustic develops 
on the source's side of the observer's lightcone. The caustic is traced by cusps in the wavefronts, as shown in the 
bottom panel of the figure. 

The equation for the time delay At (the difference i n arr ivai time between a lensed light ray and an unlensed one) 
can be obtained from our approach by solving for t in (46b) and subtracting from it the travel time in the absence of 
the lens: 



At 



y-yi 



CC0S((y5 



a 



Tl+T 



ir 



(63) 



The equation for the time delay is usually given up to second order in the angles. We need expansions of xi and Ti, 
from Eqs. (43) with r — 0, up to second-order terms: 



Vi 

Ti = — 

C 



V 2 

1 + T 



o(p 3 



xi = yi<p + 0(tp 3 ). 



(64a) 
(64b) 



These, as well as (|62|) for x, are used in (p3l). We obtain 



A T i AAa 2 

AT = T + ^D7 a > 



(65) 



which agrees with the usuai expression for the time delay. 



B. 3+1 static lens planes 

The generalization of our scheme for the thin-lens spacetimes to 3+1 dimensions is straightforward. We considcr 
a spacetime with coordinates (t,x,y, z) consisting of two Minkowski spacetimes matched together at the lens piane, 
which in this case is the timelike 3-surface z = zi and is parametrized by (t, xi, yi). A generic nuli direction is l a given 



in terms of stereographic coordinates as in Eq. (33). For our present purposes, we relabel the components of l a as 



l a = -y=(l,Px,Py,Pz) (66) 



where (p x ,p y ) label the sphere of nuli directions and p z = «/ 1 — p 2 x — p y . As in the (2+l)-case, we choose our family 
of nuli surfaces as nuli planes on the observer's side of the spacetime, and write it parametrically as 

r = ctj> + s (67) 

with f = (x,y,z), p = (p x ,p y ,p z ) and s= (s x ,s y ,s z ), with s-p = 0, representing the nuli surface at r = 0. The 
values of s label the nuli geodesics of the p nuli piane. 

Then, if ti represents the time that it takes for a lightray to reach the point fi = (xi, yi, zi) measured from the time 
that the wavefront passes the observer, we must have 

fi = ctip + s (68) 

and thus 

n = . (69) 
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It is very useful to change the geodesie labels from s to fi by means of 

s = n- (n ■ p)p. (70) 

for use on both the observer and source side of the spacetime. Eq. ([37j) becomes 

f=(cT—(frpj)p + fi. (71) 

We thus have in this picture that the wavefront, in direction p, consists of nuli geodesics parallel to each other, 
parametrized by the values of (xi,yi) where they hit the lens piane at zi. Once a ray reaches the lens piane, it is 
delayed for a time T(xi,y{), during which it just "sits" at the lens piane, before leaving in a different direction p 1 . 
Thus, in our model, we have 

f=n (72) 

and then 

f=c(T -t-Tt-T)tf + fi (73) 

where 

p> = p- a(xi,yi). (74) 

Here we have defined a(xi,yi) in the usuai manner, representing the deviation of the nuli geodesie, which depends on 
the point at which the geodesie hits the lens piane. Equations (|7l|), (|72]), and (f73j) constitute a parametric expression 
of G(t, x, y, z,p x ,py) — 0. In this approach, the time delay T(xi, yi) is assumed to be prescribed, whereas a(xi,yi) is 
determined by the envelope condition. The envelope condition in this case takes the form 

— - — - (75) 
dp x dp y 

which is carried out implicitly as in the (2+l)-case, and which yields 

dT 

c— = -a x (xi,yi), (76a) 
oxi 

dT 

c^— = -at v (xi,yi), (76b) 
àyi 

thus reproducing the relationship between the time delay and the bending angle in the more standard astrophysical 
approach to lensing. The detailed calculations, which can be omitted, are reproduced below. 
The implicit G function on the observer's side is given by the z— component of Eq. j7l| ) 

G(T,x,y,z,p Xl p y ) = z- zi - (t- n)p z = 0, (77) 

while xi and yi are defined implicitly by the (x, y)— components of (|7l|). Solving for xi and yi from these equations 
and substituting back into G yields 

G = — (-£ + xp x + yp v + zp z ) = 0. (78) 

Pz 

On the source side we implicitly define the thin-lens G function by 

G = z — zi — (t — ti — T)p' z = 0, (79) 

where xi and yi are implicitly defined by the (x, y)— components of Eq.(|75|). The envelope of G is obtained by setting 
to zero the partial derivatives of G with respect to p x and p y , holding (t, x, y, z) fixed. Applying this prescription to 
the G function before the lens piane, we have 

BG n 

-^- = => x-^z = 0, (80a) 
op x p z 

|^ = => y-Ey z = Q, ( 80 b) 
OPy Pz 
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We note that Eqs. ( Bffl imply that 

* = » = * . (81) 
P a Pz 

Solving for and p y from Eqs. ( |S0| ) and substituting into the G function, we obtain our final expression for the 
thin-lens G function on the observer's side 

PzG = -t + xp x + yp y + zp z = —t+ ^{x 2 + V 2 + z 2 ) 



= -t± y/x 2 + y 2 + z 2 = 0, (82) 

the observer's lightcone. 

On the source's side, the envelope construction requires computing the following partial derivatives: J^, Jjp-, Sp- 
aivi from the (x, y)— components of Eq. (|73|). Since, p 1 = p — ct(xi, yi), we need the derivatives of p 1 . 
We definc 



s Px _ d P'x | d Px dx i } d Px dyi 
Sp x ~ dp x dxi dp x dyi dp x ' 

and similar quantities for ali combinations of (j>' x ,p' y ,p x ,Py)- After a lengthy calculation, we obtain 



(83) 



O.r 



d P* 1 -Px§xl-Pv^-P-Px-PvPy 

d Vl n T) {-^ (l - §Lp> x - p> x p x ) - f| (gg, + p > yPx ) } 

d P* 1 - P*f*7 " P'yWt ~ PxP'x ~ PyP'y 



(84a) 
(84b) 



and similar expressions for dxi/dp y and dyi/dp y . 

The envelope construction consists in setting the p x and pj, partial derivatives (holding (t, x, y, z) fixed) of Eq. ( |79| ) 
to zero. Explicitly, we have 



= -fr - ti - T) ( + 

' V^Px^Px dp' y S Px J 



, , dxi dyi dT dxi dT dy x ' ^ 
l> ' : ' x 5p x <>.,. dxidp x iìiji dp x J ' 



-(t-ti-T) 



V^Px^Pa dp'y^Py 



, / i9a;; <9™ dT dxi dT dyi \ ,„„, , 

+pUPx^- J -+P 3 /^- + ^-^- J - + ^-/ ì , 85b 
V dp a dp y dp a óyi dp y J 

where we have again used Eq. (|8l]). In these two equations, we substitute the partial derivatives given above and 
solve for dT/dxi and dT/dyi, obtaining 

dT 

tt~ =Px ~Px = a x (xi,yi), (86a) 
oxi 

dT 

Q~-=P y -Py = ot y (xi,yi). (86b) 

Eqs. d86| ) determine the new direction of the light ray on the source's side with respect to the initial direction of 
the light ray on the observer's side (i.e., the bending angle) as the gradient of the time delay function. Note that we 
have not made use of the small angle approximation in these calculations. 
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V. MOVING LENS PLANES IN 2+1 DIMENSIONS 



As an application of our approach, we consider a generalization of our (2+l)-results to the case of a lens piane that 
is moving in the observer's trame of reference, but is otherwise unchanging. In other words, the deflector has a given 
configuration, which does not change with time, but moves rigidly along the line of sight. The idea is to obtain the 
bending angle and compare with the bending angle by the same deflector at rest. 

Again our basic picture is a family of parallel nuli geodesics at an angle ip with the normal to the lens piane forming 
a nuli piane. The rays each arrive at the lens piane at a time Ti, and are then remain in the lens piane for a length of 
time T that depends only on the point at which they hit the now moving lens piane, before exiting the lens piane at 
an angle tp' = tp — a m . 

We point out that the "detention time" T of a thin rigid lens in (slow) motion is the same as that of the same 
deflector at rest. The reason for this becomes clear if we consider the situation of a glass sheet of thickness A that is 
moving. In the rest frame of the glass, light has a speed uq and takes a time to = A/uq between entering and exiting 
the glass. The time in the laboratory frame in which the glass is moving with speed v is t' — to ~~ vA/c 2 (neglecting 
higher powers inv/c) following from the fact that the length for light to travel changes due to the motion of the glass, 
and with the fact that the speed of light also changes if the glass is moving (by the Fizeau effect). Thus t' differs 
from to only by a term proportional to the thickness of the glass, but completely unrelated to the glass properties. 
In the limit in which the thickness of the glass tends to zero while the rest-frame travel time, to is kept Constant, the 
difference vanishes. 

As before, we assume that the observer lies at the origin of coordinates Xq — (To, 0, 0). The lens piane, however, is 
moving, being described by xf = (t,Xi,yi(t)). The (affine) parameter r describing the evolution of the nuli geodesics 
of the nuli surface is such that the wavefront r = passes by the observer at f = To. The parameter r runs backwards 
in time, so that the time coordinate of the lightrays is given by t = To — r as before. We will slightly abuse our 
notation by using yi(r) to denote yi(t = T — r). If the individuai lightrays in the beam are labeled by r, with r = 
being the ray that passes through the origin, then we can give them parametrically as 

x = r coup + ct sin ip, (87a) 
y = —r sin tp + ct cos tp. (87b) 

Individuai lightrays reach the lens piane at time Ti, at which 

xi = r cos tp + cti sin tp (88a) 
Ui( T i) — —rsinp + cti cosp (88b) 

The rays remain at the lens piane for a time T{x{) and then leave in a direction tp — a m {xi), where the bending angle 
for a moving lens, a m (xi), is to be determined. At the time the rays leave the lens piane, however, the lens piane has 
moved to a point yi(ri + T), carrying the rays with it. On the other side of the lens, the source side, we have 

x = c(t — ti —T) sin(tp — a m ) + xi, (89a) 
y = c(T-T t - T) cos(<^ - a m ) + yi{n + T). (89b) 



Eqs. fl87p hold when t < t; , whcreas Eqs. ( p9[) hold for r > n+7. On the observer side (|87b| ) represents G(x, y, t, tp) = 



if r is thought of as a function of (t, x, t p) gi ven implicitly by (87a). The envelope condition dG/dp — on the observer 



side, implemented by taking d/dtp of ( S7a ), yields, as we had earlier, r = or 



xi = cti sinyj, (90a) 
Vlin) = cri cos tp. (90b) 

As earlier, we consider ( |39bj ) as the equation G(x,y,t, tp) = on the source's side by interpreting xi as a function of 

(T,x,tp): 

xi = xi(t, x, tp), (91) 

given implicitly by ( |S9a| ). The implementation of the envelope condition, dG/dtp = 0, involves several steps that must 
be explained. 

First, we have to evaluate the velocity of the lens when the light first enters it, at 77 + T and when it exits the 
lens, at ti, on the observer side. Since the lens is thin and moving slowly, we consider the two values to be the same 
(or equal to an average) and cali 
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dyi 
dr 



-vi- 



(92) 



The second point is that, though Eq. (89b) contains xi via a and T, it also contains ti which is a function of both 
xi and tp and is thus needed in the envelope condition. It is obtained by eliminating r from Eqs. (|S7l): 



cri ~ xi sintp- ìji{ti) = 0. 
Taking the tp derivative of Eq. (|9^) and using (p2[), with r — 0, yiclds 



sin </? 



<9</3 d(f (c + COS </Wj ) 

Next we need the (/5—derivative of xi, calculated from ( |S9a| ) , which, using 

dxi _ -c(r -ti-T) cos(p - a m ) da m _ ( &T 
dtp 1 - c(r - ti-T) cos(<^ - a rn ) dxi 

Finally, setting dG/dtp = from (B9q) yields 



yields 



+ 



sin tp 



dxi (c + vi cos tp) 



sm(tp - a m ). 



(93) 



(94) 



(95) 



dri dT dxi 

dtp dxi dtp 

+c(r - ti -T) sìn(tp - a m ) I 1 - 



(ceos(( y 9 - Oro) + vi) 

da m dxi 



which, using (53), (E 







and the fact that xi cos tp — sva.tpyi(ri) = (from Eqs.(pO|)) becomes, 



sin(<y9 - a m ) - (c + vi cos(p - a m )) 



(dT 

\ dxi (c + % cos tp) 



snxtp 



= 0. 



(96) 



(97) 



Now, using the rcsults of the envelope condition on the observer side where tp is a function of xi given implicitly by 
eliminating t; from Eqs. (90), we can interpret Eq.(|97|) as an equation defining the bending angle, a m , as a function 
of Xi if T(xi) is known. In the regime of small angles, keeping only linear terms in tp and a m , it reduces to 



tp - ov. 



(c + Vi) 



( dT 

\dxi ' {c + vi) 



'•P 



= 



or, the equation for a m in terms of the gradient of T: 

= - |1 + 7 



9T - M 
Ctt— = I 1 

oxi 



(98) 



(99) 



This is the relationship between the bending angle and the gravitational time delay at the lens piane when the lens 
is moving. It represents a correction of a factor of (1 + v/c) with respect to the bending angle for the same deflector 
at rest. This is in agreement with aberration ||, according to which the angles made by lightrays as observed in a 
moving frame are corrected by a factor of (1 + v/c) with respect to the angles made by the same lightrays as observed 
in the rest frame. In our case both the directions of incoming and outgoing lightrays are corrected by the same factor, 
and consequently their deviation a is corrected by the same factor. 
The lens equation 



%i + {y- Vl{n + T)) tan(y> - a m ), 



(100) 



is obtained from by eliminating t — 77 from (p9[). We are interested in the small angle regime so, neglecting higher 
order terms, we have 



xi + {y - yi)(tp - a m ), 



(101) 



where we have used the assu mptio n that T i s small (of the same order as tp. and a m ) and introduced the notation 
yi = ui(ti). From the ratio of (90a) over (90b) and for small tp, we have 
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Xl_ 

m 



(102) 



which leads, via Eq. (101), to 



x = —xi - (y — yi)a r , 
Vi 



(103) 



£ — fi 
y vi 



y-yi 
v 



(104) 



Since in fiat space the coordinate distance is equal to the angular diameter distance, we can use the substitutions 
y = D Sl % = Di and y—yi — Di s to denote the distance to the source, the distance to the lens and the distance between 
the lens piane and the source piane at the time the wavefront passes by the lens piane. With these substitutions the 
lens equation becomes 



Xì_ 

D, 



Eli 
D s 



(105) 



Thus the changes in the lens equation from lens motion amount to a factor of (1 +U;/c) in the bending angle (from 
Eq. (]9^)) and the evaluation of the instantaneous position of the lens piane at the average time that the wavefront 
reaches it. 



The time delay can be obtained by solving for r in (B9b) and subtracting the travel time of a ray in the absence of 
the lens. We have thus 



a 1 V - Vi (n + T) 

At = h n 

c cos((p — a) 



■T--y/¥T\ 



(106) 



From Eq. (105) and Eq. (90b), for small p, 



Ti 



cernii + 2^ 2 ) 



(107) 



but keeping up to second-order terms, we find the time delay equation corrected for a moving lens; 



Ar = (l + vi/c)T+^lla 2 . 



(108) 



VI. COSMOLOGICAL THIN LENSES 



In sections IV and V, we derived the implicit Fermat potentials for lensing by static and moving thin lcnses where 
the lens piane separated two regions of fiat spacetimes. In this section, we expand these results to consider lensing in 
spacetimes that (apart from the lens) are described by Friedman-Robertson-Walker (FRW) metrics. Our approach is 
to utilize the fact that ali the FRW spacetimes are conformally fiat to rescale physical quantities in the fìat lensing 
case into physical quantities for observers living in a FRW universe. We begin by noting that there is, for each 
FRW spacetime (k = ±1,0), a coordinate transformation between the "naturai" FRW coordinates (t,x,0,c/)) — x a F 
to "Minkowski" coordinates (t,x,y,z) — x a M , 

= Xm&f), (109) 

such that the metric can be written as 

dZ 2 = dr 2 - a 2 (T)dll{x,6,(j)) = n 2 (dt 2 - dx 2 - dy 2 - dz 2 ). (110) 

Here, di 2 is the metric of a Constant curvature three-space (sphere for k = +1, hyperboloid for k — —1, or fiat space 
for k = 0). In general, the conformai transformation is such that 

fi = fk(x a Hr), (111) 
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for some function Because conformai transformations preserve lightrays, the implicit Fermat potentials for lensing 
in our previous sections are implicit Fermat pot entia ls for the cosmological spacetime and can be expressed in terms 
of either coordinate system, x F or x% , via Eq. ( |109| ) . 

Often, when a symbol is used for a quantìty described in the Minkowski space and we want to dìstinguìsh it from 
the same quantity in the cosmological space we will indicate that by a twiddle, Le., A versus A. Also, often when we 
are ref ering to a quantity associated with a moving lens we will indicate that by a "fiat", e.g., the yi position of the 
moving lens becomes yi (r) = yi . 

Several subtle issues arise from the fact that, in general, a fixed space point (\, 0, <fi) ^Constant moving with the 
cosmological flow in the k — ±1 cases, has in the associated Minkowski space a coordinate velocity; thus co-moving 
lenses in the FRW coordinates are modeled by lenses that move, Le., have a coordinate velocity, in the Minkowski 
space. This, in turn, leads to the consideration of aberration affects in the bending angle and source angles which 
then influence angular-diameter distances. A second subtle issue is the relationship between the gravitational time 
delay as computed in the FRW spaces compared to those computed in the associated Minkowski space. We point out 
that the time delays are not simply conformally related; the cosmological time delay must be calculated independently 
from the Minkowski delays. This issue will be further explored later in this Section. 

We consider the universe to be two sections of a FRW spacetime appropriately matched at a lens piane. The matter 
distribution in the lens piane determines two related gravitational time delay functions: the (Minkowski) coordinate 
time delay, T, and the related cosmological proper time delay, Te, with 

T c = n,T (112) 

where fi; is the conformai factor evaluated at the lens at the time of arrivai of the ray. The relevant functions depend 
on the Minkowski coordinates in the lens piane, x;. We interpret Tc{xi) as the cosmological proper time measured by 
a person just past the lens (on the observer's side) between a lightray emitted at the source at time r s that passed 
through the lens and one emitted at the same r s that was not influenced by the lens. 

As we saw for the moving lens in Sec.V, our implicit Fermat potential, in 2+1 dimensions, is given by (|S9| ) with 
the envelope condition yielding 

dT 

a n = -(c + vi)-^— = (1 +vi/c)a, (113) 

OXi 

where Vi is the average speed of the lens during the time the lightray is delayed and a m is the bending angle observed 
by a stationary observer with a moving lens and a the bending angle for an observer comoving with the lens. Our 
lens equation (in the small angle approximation) is simply 

x = <py - (y - m)a m , (114) 

where (x, y) are the source position, yi is the position of the lens piane when the lightray leaves the lens, and ip is the 
observation angle. 

We consider how to introduce some (modified) angular-diameter distance s for lenses that are moving in a Minkowski 
spacetime. Here, we interpret the source coordinate, x, as the metric distance, A s , between the source (in the source 
piane) and the optical-axis, (i.e. x = A s ). We then convert the coordinate distances, y, and y — yi into angular- 
diameter distances, D s and D/ s , as follows. Let [3 be the unlensed angle that the source at x subtends at the observer 
and let 7 be the angle that x subtends at the lens piane (when the lightray leaves the lens piane). Then 

v = D > = 7i = ir> ( 115a ) 
p p 

y-y l = D ls = - = —. (115b) 

7 7 



Eliminating y and y — yi from ( |114| ), the lens equation becomes 

P = p-^ am . (116) 

This expression for the lens equation with a moving lens is in terms of distances and angles relative to the frame of 
a stationary observer in the Minkowski space. As we saw earlier, by the aberration (of light) equation, the bending 
angle seenby the stationary observer, a m , is related to the bending angle a seen in the frame comoving with the lens, 



by Eq. (113) 
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a(l + -), (H7) 

c 



so that Eq. (116) becomes 



= <p- — a m = tp- — a 1 + — (118) 
L>s i->s c 

with Di s the time-dependent angular-diameter distance of the lens to the source as seen by a statìonary observer. On 
the other hand, if we wish to use (as we will in a moment) the angular-diameter distance measured from the lens to 

the source in the comoving lens frame, > we have that 



D\f = -f^r = -^r, (119) 



where t he a ngle that x subtends according to an observer comoving with the lens piane. Using the aberration 



equation, Eq. (117), for the angle 7, we have 



3 * = d u ) qT*)- ( 12 °) 



The (moving) lens equation becomes 



r)( CO ) 

f3 = l p-^-a. (121) 
Note that neither d[ c °^ or Di s are directly observable quantities. The issue of whi ch to use depends on the physical 



situation that is being addressed. If we utilize the coordinate transformation, Eq. (109), that relates the Minkowski 



coordinates to the "naturai" FRW coordinates, the moving lens in the Minkowski background becomes a stationary 
lens in the conformally related cosmological spacetime. Hence, should be used in the lens equation in the case 

where our interest lies in the FRW space. 

Th e ju st completed discussion pertained to the Minkowski background; we now transform the lens equation, 



Eq. (121), to the cosmological background. We assume that the observer's worldline is given by (r, x = 0>' 
0, 4> — 0) in the standard FRW coordinates. In any of the three FRW models, such an observer is stationary at 
(t,r = 0,9 = 0,0 = 0). However, the lens and source, located along the worldlines of Constant (x, 9, <f>) will be 
"moving" in the Minkowski coordinates. (The motion of the source in the Minkowski coordinates plays no important 
role.) In the cosmological space, the metric distance between the source location and the optical axis (in the source 
piane) is given by A s = Q(x®)A s where the conformai factor is evaluated at the source location when the lightray 
leaves the source. Using 

D s = % = Q{x a s )D s (122a) 

P 

Dis = ^= n{x«)D[ C s 0) (122b) 
for the cosmological angular-diameter distances, our lens equation immediately becomes 

f3 = <p-^J?-a, (123) 

which is identical (in form) to the fìat-space lens equation with a stationary lens and to the conventional cosmological 
lens equation |l]||. 

We now turn from the cosmological lens equation to the cosmological time of arrivai equation. In Section V, we 
derived an expression for the "time delay" at the observer between the true "lensed" path and the "unlensed" path 
in a Minkowski model with a moving lens, namely, 

AT=(l + v l /c)T+^^a 2 m . (124) 
2c D s 
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Using the transformation properties of a m and Di s just described, Eq. (124) has the forni 

Ar = (l + V C )^ + ^^-a 2 j, (125) 

where a is the bending angle seen in the comoving frame and T is the Minkowski coordinate time-delay measured 
along a worldli ne ju st to the observer's side of the lens between the arrivai time for a lensed and unlensed ray. To 
transform Eq. (125) to cosmological variables, we first note that 

%- = ^ (126) 
D s D s 

where the ratio on the right side refers to cosmological angular-diameter distances. Next, we see that if A; is the 
Minkowski metric distance from the optical axis to the "image location" in the lens piane and tp is the angle this 
distance subtends at the observer we have 

5j = n^T)A L = n(a .« )5 (127) 
f 

which relates the angular-diameter distance from observer to the lens in the FRW cosmology to the conformally 
related "Minkowski" angular-diameter distance. (Here, the conformai facto r is evaluated at the location of the lens at 
the time the lightray leaves the lens.) If we multiply both sides of Eq. (|l25j ) by ilo^; where fio and fi; are respectively 
the conformai factor evaluated at the observer at the observation time and at the lens at the time the light ray leaves 
the lens, we have 

AW + */0§(T + f3w), (128) 

where we have used 

T c = SliT (129a) 
At q = fi Ar (129b) 

the latter being the cosmological proper time between reception of the two rays at the observer. 
From the relation 

(1+V<0^ = ! + (130) 
5 2/ ai 



which is derived in the appendix for small coordinate velocities, vi/c, in the associated Minkowski space, from Eq. (128), 
our final cosmological time delay equation is 

At q = {1 + z)(t c +^^cA. (131) 

Our last issue is to describe how the cosmological proper time delay, Tq is to be computed. To compute Te, and, 
by association, the Minkowksi coordinate time delay, T = VL^ l Tc, we first consider the Newtonian potential 

Uc^-frJ^-, (132) 
J \x-x'\c 

where dm' = pc(x')d 3 x' is the "observable" mass element in the cosmological space, and the norm in the denominator 
is taken in the cosmological metric. To find the cosmological gravitational time delay, one integrates the Newtonian 
potential over the path from the source to the observer, 

Tc = ~^ j Ucdlc > (133) 

where die is the distance element along the path in the cosmology. In the standard thin-lens approximation, it is 
assumed that pc{%') is non-zero only in a thin region perpendicular to the optical axis. Thus, when computing the 
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time delay from the integrai in Eq. ( |L33| ) , one need only consider the region along the trajectory that is very dose to 
the lens piane. First, we point out that under conformai rescaling, since dm' is conformally invariant, that 

/dm' f dm' „ , , 

,^ ~ - / „-p ^| = n,~ l U M 134 

\x-x'\c J \li\x-x'\ M 

where \x — x'\m = is the metric distance between two points in the lens piane taken in the conformally related 
Minkowski metric. Because Uc only has support in the vicinity of the lens and die = £Mm, where dlu is the 
distance element in a Minkowski spacetime, we have 

T c = ^ J U c dl c = Z J ì U M dl M = T M . (135) 

There is an apparent conundrum here. Te is the physical proper-time delay which appears to be the same as the 
Minkowski space proper-time delay for the same physical situation, Le., arising from potentials 

U C = - f rJ^- (136a) 
J \x-x'\c 

U M = - f :J m ' (136b) 
J \x-x'\m 

with the same dm' in both cases. The conundrum is resolved by noting that though the dm' are the same the 
density functions are different. 

We compute the mass densities by introducing physical coordinates at the lens (locai Lorentzian coordinates), so 
that the cosmological metric at the lens takes the form 

ds 2 = (dt 2 c - dx 2 c - dy 2 c - dz 2 c ) = ttf (dt 2 - dx 2 - dy 2 - dz 2 ) (137) 

with 

x a c = Q lX a . (138) 

Using 

dm' = p c (x^)d 3 x c = p c {ttix a )tfd 3 x = p M (x a )d 3 x (139) 

we have the relationship between the physical density, pc, and the equìvalent fictìtious Minkowski space density 

PM(x a ), 

pM(x a ) = pc{ttix a )nf. (140) 



Returning to the issue of the cosmological lens equation, Eq. (|123| ), and associated bending angle, Eq. (|113j ), we se e 
t hat the present results concerning the cosmological time-delay agree with the earlier results, since from Eqs. (129a), 
(138) and (113), we have 

dT dT c 
dxi Uidxi dxci 

the cosmological bending angle. 

Ali our results obtained from the envelope construction now agree with the standard results appearing in the 
published literature. In additimi, it is clear that the method can be used to extend the results to the case of a lens 
piane that has a peculiar motion above the Hubble flow by letting the velocity of the lens piane in the underlaying 
fiat space be arbitrary. 

As a final comment here, we point out that this method of obtaining the FRW lens equation and time delay equation 
by conformai rescaling of an associated Minkowski space can be extended, largely unchanged, to any spacetime that 
is conformally related to Minkowski space. It is not clear whether this observation has any physical application. 
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VII. CONCLUDING REMARKS AND OUTLOOK 



We have introduced a novel approach to gravitational lensing built on a variational principle that is analogous to, 
but distinctly different from, the conventional version of Fermat's principle. 

Our approach is to find the envelope of an appropriate family of nuli surfaces containing the observer's location. 
We have shown how this approach reproduces the astrophysical scenario of static thin lenses as an illustration. More 
interestingly, though, we have been able to use this approach to obtain the correction to first order in v/c that 
the motion of the deflector imposes on the bending angle in the approximation of thin lenses. This calculation can 
alternatively be done by direct integration of the nuli geodesics in time-dependent spacetime perturbations off fiat 
space. In this respect, our result agrees with such a direct calculation of the bending angle by moving thin lenses 
carried out by Pyne and Birkinshaw whereas it differs from others jlO| by an overall sign and a factor of 2. Work is 
in progress to clarify this difference and will be reported elsewhere. Notice that a redshift of 0.001 might conceivably 
bring this effect into the observable regime in the future. 

Because in the thin-lens case in cosmology our approach is based entirely on the under lying conformally fiat space, 
we also have an alternative derivation of the lens equation in cosmology, which does not seem to have been exploited 
in the literature so far. We feel our derivation clarifies some points that remain obscure in the presentations of the 
lens mapping in cosmology that appear in |l|] and ||. On the other hand, the reader interested in a very complete 
derivation of the lens mapping in cosmology by integration of the nuli geodesics will definitely benefit from the excellent 
article by Pyne and Birkinshaw JÌT| . 

Throughout this work, we have assumed that a metric is given that represents the structure of the deflector or lens. 
In this sense, we have kept ourselves within the kinematics of the implicit Fermat potential. We have not addressed 
the issue of the dynamics of the Fermat potential, namely: how the implicit Fermat potential is directly affected by 
the structure of the deflectors. The resolution of this important question in principle involves two steps: to solve 
the Einstein equations for the metric, and then to use the metric as a given source into the eikonal equation for the 
implicit Fermat potential. By contrast, in the standard thin-lens scenario, the Fermat potential is obtained directly 
from the surface mass distribution by means of a 2-dimensional equation of the Poisson type. It is reasonable to 
ask whether an analogous scheme to obtain field equations for the implicit Fermat potential directly in terms of the 
structure of the deflector exists in cases other than the thin-lens scenario. This question will be discussed elsewhere. 
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APPENDIX: 



In the main text we quoted the approximate equation, 

+ (Al) 

which relates the ratio of conformai factors fi to the cosmological scale factor a, at two different points, the observer 
and the lens, that are connected by a nuli geodesie. The quantity v\ is the velocity of the lens in the associated 
Minkowski coordinates evaluated as the ray leaves the lens. 

For the case of k = this result is trivially true since then a = il and vi = 0. For ease of presentation we will only 
give the proof for the k = — 1 spacetime, but the calculations are exactly paralled in the case of k = +1 and the result 
is exactly the same. We start by recalling that we have the metric in the two forms 

(tì 2 = a 2 (dT] 2 - d X 2 - sinh 2 X (d0 2 + sin 2 0dtp 2 )) (A2a) 
= n 2 (c 2 dt 2 - dr 2 - r 2 {d0 2 + sin 2 Odip 2 ) ) (A2b) 



et = — '— — (A3a) 

cosh rj + cosh \ 

SÌnhx (A3b) 



with Eq. (JL09|) explicitly given by 



cosh 77 + cosh x 
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and 



r> 2 = 



Aa z 



(1 - (et + r) 2 )(l - (et - r) 2 )' 



or, from Eq. ( |A3|) , 



n 



coshx + coshfy. 



(A4) 



(A5) 



We are interested in a lens at some fixed value of 9, <p and x — Xi > sending a light signal to the observer located at 
X = 0. The lens and arrivai times are r\i and i] , respectively, and we have 



m-vi = xi- 

Using Eqs. ([A|) and ( |A6| ) we can construct 

ili _ tt coshx; + cosh(x; - rjo) 



or, for small Xh 



By calculating the velocity from 



we have 



ai a 

fi; ^ ilo 

ai ao 



1 + cosh rjo 

Xi sinh rjo 
1 + cosh rjo 



dr_ _ %j\x 
* " S Ix 



sinh % sinh 7] 
1 + cosh % cosh ?7 



(A6) 
(A7) 

(A8) 

(A9) 
(AIO) 



We thus see that a FWR lens piane approaches the observer in the conformally related Minkowski space. Evaluating 
at the lens time we have 
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sinhx/ sinh(x; - m) 



and for small values of xi 



c 1 + cosh xi cosh(xz - ?/o) 
Xi sinh rjo 



!1 

c 



1 + cosh rio 



(Ali) 



(A12) 



Using Eq. (A12) with (A8), we have our result, Eq. (Al). 
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FIG. 1. Parametric construction of the nuli surfaces. On the observer side, the surface is a nuli piane ruled by parallel 
lightrays making an angle ip with the y— axis. The lightrays are labeled by their distance r away from the observer, and are 
parametrized by the distance cr away from the wavefront that passes by the observer. On the other side of the lens, the 
lightrays bend through an unspecifìed angle a that depends on the position at which the lightrays reach the lens piane. 
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FIG. 2. Envelope construction for the case of a bending angle independent of the location on the lens line xi. The top panel 
shows a spacetime picture of one member of the 1-parameter family of nuli surfaces that is used to construct the envelope. 
The center panel is a space picture of the nuli surface, where both the lightrays and the associated wavefronts are shown. The 
bottoni panel shows a spacetime picture of the envelope of the 1-parameter family. The family is labeled by the direction of 
the piane wave that hits the lens line. 
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FIG. 3. Envelope construction for the case of a time delay T{xì) = Tq — hi(l + (xi/xo) 2 )- The top panel shows a space 
projection of one member of the family of nuli surfaces. The nuli surfaces are labeled by the direction of the piane that hits the 
lens line. The center panel contains a space picture of the envelope showing the lightrays. Notice the caustic line that develops 
after the lightrays pass through the lens line. The bottoni panel shows a greatly magnified view of a few of the wavefronts 
associated with the lightrays of the center panel. Notice the cusp in each wavefront. These wavefronts occur around the value 
20 in the verical axis of the center panel. 
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